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Abstract 

The sensitivity of an ice radiowave detector to the anticipated high energy neutrino fluxes 
is calculated on the basis of a detailed threshold analysis and computation of the shower 
production rate. We show that diffuse neutrinos from Active Galactic Nuclei could be 
detectable in a radio detector of 1 km 2 area established in Central Antarctica. 
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About 30 years ago Askar'yan proposed a new method for the detection of high energy particles 
ON by means of the Cherenkov coherent radiowave emission from the negative charge excess of 
electromagnetic showers generated in air or dense media []]]. The charge imbalance of a shower 
is created by the Compton scattering of shower photons on atomic electrons, the annihilation 
Ph! of shower positrons in flight and the knock-on process. The percentage of negative charge 
excess amounts to ~ 20% at the shower maximum. The resulting Cherenkov emission by 
excess electrons is coherent at wavelengths larger than the shower lateral dimension, i.e. in the 
radiowave region. In spite of the very low frequencies compared with visible light, this emission 
should be observable for sufficiently high primary energy because the radiated power scales 
with the square of the shower size. 

The radio pulses have been indeed successfully observed from air showers (EAS) in coinci- 
dence with particle arrays, but some their properties (polarization, pulse shape, south-north 
asymmetry of pulse amplitudes, etc.) indicate that the radio emission has a different dominant 
mechanism of generation in air. Nevertheless, by theoretical reasons the shower charge excess 
should be the principal source of the coherent radio emission in dense media, and after a while 
interest in the idea of Askar'yan was renewed by the suggestion to detect high energy (HE) 
neutrinos in cold antarctic ice, which has very low radiowave absorption at temperatures below 
— 50°C 0. It was argued that a radio antenna array placed on the glacier surface in Central 
Antarctica could provide an effective target volume of the order of 10 9 — 10 10 m 3 for cosmic 
neutrinos with energies above 100 TeV 0— 

In this paper we make a detailed threshold analysis for such a detector, taking account 
of radiowave absorption in ice. We calculate as well a production rate of showers initiated 
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by upward-going HE neutrinos in the upper layer of the Earth. That allows us to estimate an 
expected neutrino event rate in a radiowave detector from the different anticipated HE neutrino 
fluxes. 



2 Event rate in a radio detector 

We consider a radiowave antarctic neutrino detector as a number of downward-looking antennae 
disposed on several dozen meters triangular grid enclosing a glacier area of the order of 1 km 2 . 
The antennae sample Cherenkov radio pulses from ice, which would provide well-defined conic- 
type images on the grid (the Cherenkov angle is equal to 56° in ice) @H(f]- 

In general, the radio detector event rate induced by the neutrino flux is given by 

N v = J dE J dQV eff (E ,e,<P) Pu (E ,e,<p). (1) 

Here p u (E , 6, 0) is the electromagnetic shower production rate per unit volume for a given 
shower energy E D and direction (9,<f>), V e ff(E ,9,4>) is the effective detection volume. 

To make an initial estimation of the event rate we neglect the angular dependence of V e ff 
and use the next approximation: 

Veff = S[R max (E ) — Rmin), (2) 

where S is the enclosed grid area, R max (E ) is the maximum detection depth for a given shower 
energy E Q , R m i n is the minimum detection depth, which depends mainly on the grid spacing. 
In this case Exp. (1) is simplified: 

N v = J dE V eff (E )P u (E ), (3) 

>E 

nit n 

where 

p v {e ) = J dn Pv (E ,e,<f>) 

and the minimum detection energy E min corresponds to the minimum detection depth R m in- 



3 Threshold analysis 

A numerical real time computation of the radiowave emission from electromagnetic showers 
developed in the totally transparent ice || gives the following parametrization of the electric 
field spectrum at the Cherenkov angle: 

Here v is the frequency, R is the distance from the shower, E a is the incident electron (photon) 
energy, v Q = 500 MHz (the result of 0] was divided by 2 to define the Fourier transform as 
E{lo) = J dtE(t)exp(iujt)). 
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Radiowave attenuation in the real ice strongly depends on the wavelength, as well as on 
the ice temperature. Therefore, we have used the original data on radiowave absorption in ice 
H together with the results of temperature profile measurements in a super deep bore hole at 
the Vostok Antarctic Station |7| to calculate the total attenuation of a shower radio pulse after 
vertical propagation from a given depth to the ice surface (Fig. 1). 

To obtain the threshold energy for the one-channel radiowave detection of electromagnetic 
showers we need to consider in some detail a process of radio pulse transformations by a receiving 
antenna and its preamplifier (active filter). The relation between output antenna voltage V(u) 
and incident electric field E(uj) is given by 

V(uj) = R a (uj)-E(u J ), (5) 

where Ra{w) is the reception transfer function of antenna. As it is recognized, the so-called 
TEM horn is the most promising broad-band antenna for impulsive field measurements (for 
instance, see Ref. ||). Its reception transfer function has approximately constant magnitude 
and linear phase dependence over the frequency range from a hundred MHz to several GHz. 
As one can see (Fig. 1), the shower radio pulse spectrum has the same band of wavelengths. 
Therefore, the TEM horn will produce an output voltage that is a high fidelity replica of the 
shower electric field in the time domain. For example, the lxl m 2 arcsine TEM horn specially 
designed for the neutrino radiowave experiments has Ra — 0.14 V/(V/m) for a normal incident 
field in air. The reception pattern of this antenna is rather broad with the half-amplitude 
beamwidths of about 90° in both E and H planes for an incident electromagnetic pulse of 1 ns 
duration 0]. 

The conventional definition of a signal-to-noise ratio at the filter output is || 

S peak instantaneous output signal power 
N output noise power 

The maximum of (6) occurs when the filter transfer function is proportional to the complex 
conjugation of the input filter voltage V*(u) ("matched filter"). The maximum value is 

[S/N] max = -J f^M| 2 , (7) 

— oo 

where N a is the one-sided white noise power density at the filter input. From (5) and (7) we 
obtain for the normal incident shower radio pulse received by a TEM horn in ice: 

[S/N] max = —eR\ J ^-\E(u,R,6 d )\\ (8) 

— oo 

where e is the relative permeability of ice (Ra rises by a factor of \fe in medium). If the antenna 
impedance is equal to the load (filter) one, 

N = kT N Z L . (9) 

Here k = 1.381023 J/K is the Boltzmann constant, T/v is the noise temperature and Zl is 
the load impedance. Using the parametrization (4) and taking into account the radiowave 
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absorption in ice for vertical pulse propagation from 100 m depth, we find for Zl = 50 Ohm 
and \fe = 1.8 (ice refraction coefficient): 

[S/N] max ^ o.i^-pg- [ M ^) . (10) 

For the signal-to-noise ratio of unity, Tjy = 300 K (according to the antarctic noise measure- 



ments [ID]) and Ra = 0.14 V/(V/m) the shower threshold energy E t h is approximately equal to 
55 TeV. This is a factor of 7 lower than the result obtained by Zas, Halzen and Stanev ||. They 
performed the threshold estimation for a half-wave dipole antenna adapted for the narrow-band 
receiving technique of EAS radio detection experiments. Hence, the use of a broad-band TEM 
horn together with a matched filter result in the significant decrease of the threshold energy. 

Fig. 2 shows the calculated dependence of the relative threshold energy on the shower pro- 
duction depth. The related results of the maximum detection depth R max can be parameterized 
by the expression: 

RmaAEo) ^ 615 ln(l + E o /5.8E 100 ) (m), (11) 
where E wo is the threshold energy for the depth of 100 m. 



4 Production rate of neutrino-induced showers 

Electromagnetic showers are directly produced only by the charged current interactions of 
electron neutrinos: 

u e + N — > e + . . . . 

\ (12) 
e.m. shower 

But hadronic showers, which are generated by all the types of neutrino interactions 

v + N — > I + hadrons , 

\ (12') 
shower 

initiate electromagnetic subshowers too through decay of neutral pions and eta particles. Ac- 
cording to the accelerator calorimetric experiments and Monte Carlo simulations, the contri- 
bution of electromagnetic subshowers to the total energy deposit of hadronic shower increases 



with the initial hadron energy and amounts to ~ 90% at E Q > 10 TeV [TT[. As the mean 



inelasticity in the reaction (12') is ~ 0.3 for E v > 100 Tev [12], one should expect that the 
power in the radio signal from the neutrino-induced hadronic shower will be only slightly less 
than the power from the pure electromagnetic shower of the same energy. In the paper we 
neglect this difference, though the question deserves more detailed study. 

The shower production rate P J/ (£', 7) for upward-going neutrinos is defined as 

P v (E ,i) = n N o G™ {E 0i i)*{E ,i). (13) 
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Here G t ° t {E 0) r j) is the total shower production function, $(i£ , 7) is the differential power-law 
energy spectrum of incident neutrino flux with the integral spectral index 7, is the nucleon 
number per unit volume; a — 1.1 x 10~ 34 cm 2 is the standard cross section. 

The total shower production functions G t ° t {E 0) 7) for electron and muon neutrino fluxes are 
given by 



G%{E„, 7) = ^^n eff (E , 7) + G?(E , 7) (14) 

G%(E 0} 7) = Gt c (E 0} 7) + G?(E , 7) (14') 

respectively with the definition of G 3 v (E ,j) (j = cc for the charged current interaction and 
j = nc for the neutral one) as 

G{{E 0ll ) = J dyy^ eff (Ejy n ) daJ ^jy> y) (15) 

In these expressions a^{E) is the neutrino-nucleon cross section, y = (E v —Ei)/E u is the relative 
energy loss of neutrino in the laboratory frame and VL e ff{E^) is the effective solid angle for 
upward-going neutrino flux, which is defined through the coefficient of neutrino absorption by 
the Earth k(E,9,^) as in (T3| 



1 

tt eff (E,7) = 2tt J dcos6k{E,6,j). (16) 


Let us notice that the shower production functions Gi{E ,^) are actually the generalization of 
the neutrino spectrum-weighted moments 

(17) 

J a Q dy 

which have been originally calculated in the paper [|TJ|] , for the case of neutrino absorption. We 
performed the computations of G^^Eo, 7) by use of the EHQL |15| and DLA |1| parametriza- 
tions of the quark distribution functions, as it has been proposed for the neutrino cross section 
evaluation in |T7|, [12] , and PREM model of the Earth structure [[L8| . The results for electron 
and muon neutrino fluxes with 7 = 2.1 are given in Fig. 3 and Fig. 4 respectively. 

All the functions GKEq,^) vary rather slowly with 7. The behavior of the v e shower 
production function around 6.4 PeV (Fig. 4) is strongly affected by the Glashow resonance 
reaction 

v e + e" -> W~ -> all (19) 
with a very high cross section [ |T9"| , fZH| ]. The Earth is almost opaque for such neutrinos. 
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5 Detection of AGN neutrinos 



Active Galactic Nuclei (AGN) have long ago been considered as possible cosmic objects for 
high energy neutrino production |21, ^2|. Due to the recent development of the idea in some 
important ways by Stecker et al. p3|, E23J] and other authors 



24|] and other authors [P5[|- [P?|1 , nowadays AGN are 
recognized as most promising sources of HE neutrinos. Therefore, we calculate the expected 
event rate in a 1 km 2 radiowave detector for the different existing models of AGN neutrino 
production, assuming E wo = 55 TeV and R m in = 100 m (for comparative discussion of the 
models see Ref. The results are given in Table 1. The expected background from 

atmospheric neutrinos 



30] will be ~ 20 events per year. 



Table 1. Expected event rate for the different models of AGN neutrino production 



Diffuse neutrino fluxes 


Expected event rates (S/N = 1), 
year 1 km 2 


from AGN by 


Stecker et al., 1992 


110 


Szabo and Protheroe, 1992 (max) 


2800 


Szabo and Protheroe, 1992 (min) 


730 


Biermann, 1992 


110 


Sikora and Begelman, 1992 


100 



These figures can be compared with the predicted muon event rate in a 10 4 m 2 optical 
underwater neutrino telescope, such as DUMAND II, NESTOR or NT-200, which are under 
construction now. According to |29| at muon energies above 1 TeV the rate will be for Ref. [^4[ 
~ 30 per year and for Ref. p5f from 160 to 800, decreasing by several times at 10 TeV. Ultra 



high energy neutrinos from AGN can be detected by DUMAND II due to the Glashow resonance 
reaction with the rate for Ref. |2~1| about several dozen events per year [SI]. 

Fig. 5 shows a plot of rate in a radiowave detector versus signal-to-noise ratio. As one can 
see, even at high detection thresholds the diffuse AGN flux remains measurable by a radiowave 
detector for all the models, save the Sikora and Begelman one. 

We have estimated the event rate in a radiowave detector only for upward-going neutrino 
fluxes, though inclined neutrinos from the upper hemisphere would give rise to the event rate 
too, if some part of initiated Cherenkov radio emission has a direction toward the ice surface. 
However, these calculations are very sensitive to the particular form of the angular dependence 
of the effective detection volume and require a full Monte Carlo simulation of the detector 
response. 

Besides of the surface deployment of a radio array, it was proposed also to put antennae 
deep under the ice in bore holes to detect neutrinos from the upper hemisphere more effectively 
||32|| . Because of the limiting hole size (20" diameter), only small conical antennae are suitable 
in this case and, therefore, the detection threshold should be several times higher. However, a 
downward-going HE neutrino flux does not suffer from the Earth shadowing effect and this can 
partially compensate the event rate decrease, especially due to the Glashow resonance reaction. 
Some type of a combined detector is possible also. 
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6 Conclusions 



In this paper we have shown that a 1 km 2 radiowave neutrino detector established in Central 
Antarctica (at the South Pole or Vostok Station) should be sensitive to the predicted diffuse 
fluxes of AGN neutrinos at energies above a hundred TeV. For some production models AGN 
neutrinos would be effectively detected in the broad energy region up to 1 PeV, that gives, in 
principle, a possibility to determine the spectrum shape. 

The performed calculations of the shower production rate by upward-going HE neutrino 
flux are relevant also for the other HE neutrino detection experiments. 
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Figure 1: Radio pulse spectrum after vertical propagation from a given depth (the maximum 
value is equal to 43.5 nV/MHz/TeV). 
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Figure 2: Relative threshold dependence on depth (£100 is taken to be unity). 
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Figure 3: Shower production functions for muon neutrino (solid line) and antineutrino (dashed 
line). Here CC denotes the charged current and NC the neutral current. 




Figure 4: Shower production functions for electron neutrino (solid line) and antineutrino 
(dashed line). 



10 



D 




1 23456789 10 

S/N ratio 



Figure 5: Event rate versus signal-to-noise ratio (S/N = 1 for E wo = 55 TeV) for the different 
models of AGN neutrino production (SDSS, Stecker et al. (1992); SPmax and SPmin, Szabo 
and Protheroe (1992); Bi, Biermann (1992); SB, Sikora and Begelman (1992)). 
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